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We consider the sector of Horndeski’s gravity characterized by the coupling between the kinetic scalar 
field term and the Einstein tensor. We numerically construct neutron star configurations where the 
external geometry is ideirtical to the Schwarzschild metric but the iirterior structure is considerably 
different from standard general relativity. We constrain the only parameter of this model from 
the requirement that compact configurations exist, and we argue that solutions less compact than 
neutron stars, such as white dwarfs, are also supported. Therefore, our model provides an explicit 
modification of general relativity that is astrophysically viable and does not conflict with Solar 
System tests. 


Introduction. — Compact objects, such as neutron 
stars are likely to unveil soon an unexplored corner of 
gravity. Indeed, the upcoming gravitational wave obser¬ 
vation facilities are expected to probe the strong field 
and strong velocity regime of gravity typical of these ex¬ 
treme astrophysical objects [ll l3- To date, the theory 
of general relativity (GR) has passed all tests with flying 
colors, therefore it is very difficult to come up with ex¬ 
tensions that do not violate any of these tests. On the 
other hand, modifications of GR are often invoked to cope 
with quantum effects at high energy and large curvature, 
in particular in the presence of spacetime singularities. 

From a phenomenological point of view there are many 
ways to modify gravity and astrOThysical data proved 
useful to assess their viability [3, |^. Historically, one of 
the first extension of GR consists in the addition of a 
new dynamical degree of freedom in the form of a scalar 
field minimally coupled to gravity . After this pioneer¬ 
ing work, the so-called tensor-scalar theories of gravity 
developed enormously 

In this paper, we consider a particular extension of GR 
that was originally discovered by Horndeski in the early 
seventies @ consisting in the most general scalar-tensor 
theory in four dimensions with at most second order dif¬ 
ferential equations of motion. This feature is particularly 
desirable since, while several theories with higher deriva¬ 
tive terms exist (see e.g. 0), they usually suffer from 
ghost instabilities or fail to be hyperbolic 0 (see also 0 
for alternative models of gravity without ghosts). 

Recently, Horndeski’s gravity was rediscovered and 
called Galileon gravity, the name originating from the in¬ 
variance of the equations of motion under arbitrary shifts 
of the scalar held in Minkowski spacetime 0- In fact. 


the covariantized version of this theory was proven to be 
equivalent to Horndeski’s theory in [lll |. An important 
subclass of Horndeski’s gravity shows self-tuning proper¬ 
ties offering a natural explanation for the cosmological 
constant. This subclass, named “Fab Four” in [l0, con¬ 
sists in the combination of four main actions, each with 
different properties that are in turn relevant for cosmol¬ 
ogy or astrophysics. 

One particular combination of these actions, named 
“John -|- George”, has the form 

S= j ,/^d‘^x[K{R-2K) 

- ^+ Sm , ( 1 ) 

where ^ is a real scalar held, is the Einstein tensor, 
K = l/(167rG), being G is the Newton’s constant, a,r] 
are real parameters, A is the cosmological constant, and 
Sm is the action for ordinary matter helds, supposed to 
be minimally coupled to gravity in the usual way. As 
we will briehy review below, this model has exact an¬ 
alytical black hole solutions, which are asymptotically 
anti-de Sitter and have a non-trivial scalar held distri¬ 
bution. Therefore, in order to be a viable alternative 
to GR, the deviations from asymptotic hatness needs to 
be small enough to pass at least Solar System tests. A 
hrst post-Newtonian analysis of this model (with vanish¬ 
ing cosmological constant) was presented in 0 , where 
very stringent constraints on the parameters of the the¬ 
ory were found. 

Here, we are interested in spherically symmetric solu¬ 
tions of this theory in the presence of matter, focussing in 
particular on neutron star models. Our goal is to investi- 
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gate the effect of the derivative coupling on the structure 
of neutron stars and to constrain the parameter space. 

While the lack of asymptotic flatness seems to make 
these models physically unrealistic, it turns out that 
there exist a region of the parameter space such that the 
exterior metric is exactly the same as the Schwarzschild 
one because the scalar field does not backreact on the ge¬ 
ometry. On the opposite, inside the star the scalar field 
behaves as a sort of “internal hair” that interacts with 
the metric and modifies pressure and matter density dis¬ 
tributions considerably with respect to GR. The size of 
these effects depends on one parameter and the interest¬ 
ing aspect is that the deviations from GR persist also 
when this parameter vanishes. Therefore, at the level of 
the equations of motion, there is no smooth limit to GR. 

We begin by reviewing the black hole solutions of dD) 
and their properties. We then consider the modifications 
induced by the presence of a matter fluid and study the 
configuration with external Schwarzschild metric. Next, 
we solve the equations numerically in the range of mass 
and density typical of neutron stars and compare the 
results with GR. 

Black hole solutions. — The equations of motion de¬ 
riving from the action (El) read 


, 

= 0 , 


( 2 ) 

(3) 

(4) 


where 
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( 6 ) 
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and where u is the unit 4-velocity of a perfect fluid with 
pressure P and density p. 

One of the first analytical black hole solutions was 
found in 1^, by imposing the diagonal metric 


, where, from now on, the prime denotes the deriva¬ 
tive with respect to r. In addition, to avoid a diverging 
current norm on the horizon, it is necessary to set J” = 0 
[l5| . In this configuration, becomes negative outside 
the event horizon, but this seems not to introduce any 
thermodynamical instability [l3 |. This model was fur¬ 
ther developed in a series of papers, see e.g. |16l424| . In 


particular, it was shown in m that the scalar field can be 
real outside the event horizon provided one adds a nega¬ 
tive cosmological constant A in the action. In this case, (j) 
is finite at the horizon but (j)' diverges. The stability and 
the quasi-normal mode spectrum of these asymptotically 
AdS solutions were investigated in 22. [^|. 

A more general solution was proposed in [13 , where 
the scalar field has the form r) = Qt + F{r), for some 
function F and constant Q. In this case, the equations 
do not depend explicitly on time and J” = 0 is required 
to satisfy the (tr) metric equations. Note that this sector 
of Horndeski’s gravity is equivalent to a f{R,T,TpvR^'') 
model, where T^i, is the scalar field stress tensor [ 2 ^ . 

It is interesting to understand in which limit one re¬ 
covers GR. In the case Q = 0, this limit is A = 0 and 
ry —>■ 00 , revealing that ry is a non-perturbative parameter 
if a is fixed 1^. For Q ^ Q, the limit A = a = 0 leads 
to a non-trivial scalar field with a Schwarzschild metric. 
This happens because the scalar field satisfies Flp^, = 0 so 
it doesn’t backreact on the metric. Such solutions were 
dubbed “stealth” configurations in [17 1. 

The crucial point is that these considerations hold true 
only in vacuum. In fact, inside matter H^i, ^ 0 and the 
metric is affected, also when A = a = 0, resulting a 
modified theory of gravity whose deviations from GR are 
well hidden inside the star, so that they do not conflict 
with Solar System tests. In the following we focus on this 
case. 

Modified Tolman-Oppenheimer-Volkoff equations. — 
According to the discussion above, we set a = A = 0 and 
look for neutron star-like solutions assuming the usual 
equation for the matter fluid 


P' + iP + p)^=0- (9) 

The scalar current components are 

The (tt) and (rr) components of eq. ([2]) yield other two 
equations and we have 


= -b{r)df^ + -I- r'^dn'^ , (8) 

f{r) 

with = 0 and A = 0. Thanks to the shift symme¬ 
try ^ const, the equation of motion for the scalar 
field turns into the current conservation law ([31). This 
symmetry also implies that the solution depends only on 


rfb' = il-f)b, Af = -B, (10) 


where 


A = rbf ^{Pr^ + 4 :k) — 3r]Q^r , (11) 

5 = 3(1 - f)r,Q^ 

+ bf~^ [6r^/P -b (1 -b f)r‘^p - 4k(1 - /)]. 
















3 



FIG. 1. Range of QqT] where F' is real and where compact 
solutions exist. Here, Pc/Pc is dimensionless. The red (blue) 
shaded area corresponds to F'^ > 0 for rj < 0 {rj > 0). Typical 
solar mass neutron stars lie inside the black rectangle. 


Finally, the scalar field equation 


TlfbF'^ = (1 - f)7jQ^ + bPr^ , 


( 12 ) 


closes the system. By expanding around r = 0 (with the 
regularity conditions 6(0) = 6o > 0, 6'(0) = 0, /(O) = 1 
and P{0) = Pc), we find that 


’12 _ ^2 f ^ _ 2Qo(3Pe + Pc) \ 

^ V ^ 3(3Qg?7 - 4 k) J 

,^p {Pc + Pe)(3Pe + Pc) 2 
6 {iQlrj - 4k) 


+ Oir^), 
+ 0{r^), 


(13) 


where we rescaled Qq = /bo. We note that SQ^r] < 4k 
corresponds to the condition P"(0) < 0. In fact, the lat¬ 
ter is the only physically acceptable condition. One can 
show that, if this condition is violated, P is a monotonic 
growing function of r so compact configurations cannot 
exist. This constraint, however, is not sufficient to guar¬ 
antee that F'{r) is real inside the star when rj is negative. 
In fact, in this case one needs also 


1 


47r 



< Qo\v\ < 


127r 


(14) 


All these conditions are summarized in Fig. [T] 

For the numerical computations we choose the poly¬ 
trope equation of state P = where pB is the 

baryonic mass density, n is the polytropic index and 
iF is a constant. Integrating the first law of thermo¬ 
dynamics leads to p = P -I- (P/Ar)"/^"’+^), and we set 
K = I23Mq, n = 2 since this models leads to compact 
oWects with accepted mass and radius of neutron stars 
[26|. We leave a systematic study of tabulated equation 
of states (as well as cases with non-vanishing a and A) 
for future work. 

Equipped with the equation of state, we numerically 
solve following this strategy: we choose P(0) = 


Pc for a set of values for Q and integrate the system as 
a Cauchy problem up to the star’s surface at r = r*, 
determined by P{r*) = 0. Next, we match the interior 
solution to the Schwarzschild external solution and to a 
constant scalar field. The boundary conditions are 6(0) = 
1, 6'(0) = 0, P(0) = Pc- Finally, we set 77 = ±1 without 
loss of generality since the system depends on the product 
QIp only. 

The total mass M is determined by solving b{r^,) = 
6 oo(l — 2 M/r*), 6 '(r*) = 2 M 6 oo/ 7 'J, where M is the mass 
and 600 is a constant. We redefine the time coordinate 
according to tp = ty/b^, so it matches with the time 
measured by a distant observer in the flat asymptotic 
region. As a consequence, since the scalar field is lin¬ 
ear in time, we must also rescale Qoo = Q/V^oo , which 
is the value measured by the same distant observer. In 
summary, the input parameters are the central pressure 
and the bare value of Q, from which we numerically com¬ 
pute Qoo and the mass-radius relation. We work in units 
where G = c = 1. 

Note that the typical values of central density and 
pressure for white dwarfs ( pc ~ 10 ® g.cm“® and Pc ~ 
2.5 10^^ Dyne.cm”^, obtained from realistic equation of 
state tables describing nuclear matter in the low density 
regime) lead to larger values of pdPc than for neutron 
stars. Since they are within the range of parameters lead¬ 
ing to compact configurations with a real scalar field, we 
infer that white dwarfs also exists in our theory. 

Results. — One of our main results is that, for 77 < 0 
(77 > 0 ), the model allows for bigger (smaller) and more 
(less) massive neutron stars than in GR. This is illus¬ 
trated in Fig. where we show the mass-radius relation 
for different values of Qoo and 77 = ±1. Note that a sim¬ 
ilar effect was found for anisotropic stars in 23 . From 
Fig. [3] we see that the range of allowed central pressures 
shrinks with increasing Qoo. For a given value of Qoo, 
the maximum central pressure is determined by the con¬ 
straint pQlc < 4 k /3 and its value is marked by the black 
dots in Fig. [2l which mark a sharp endpoint for most 
curves. 

We remark that our solutions are not smoothly con¬ 
nected to GR. The limit Qoo = 0 does not lead to gen¬ 
eral relativity, neither does the 77 —>■ 00 limit. This was 
expected from the discussion in the last section. From 
the action it would appear that the natural GR limit is 
77 —>■ 0 , but from eq. m, we see that the scalar field 
diverges in this case. 

It is interesting to look at the scalar field inside the 
star. As mentioned above, for Q = 0, P' vanishes ex¬ 
actly outside the star. On the opposite, for Q ^ 0, the 
field smoothy connects, through the star’s surface, to the 
external field that decays at infinity, without affecting 
the metric. This is shown in Fig.[4l where we display the 
plot of F'^(r). We note that, as anticipated, the field is 
always real for 77 > 0 , while it can be imaginary for 77 < 0. 
In both cases, however, the field is real for r* < r < 00 . 
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FIG. 2. The mass-radius relation for various values of Qoo 
and »7 = -fl, —1. The thick black curve is the GR prediction 
while the dashed one has Qao ~ 0. The black dots are the 
points where the solutions cease to exist. Note that, for rj > 0, 
the curves with Qao ~ 0.016 and Qoo = 0.032 do not reach 
such a point in the chosen range of central pressures. 



FIG. 3. Lines of constant Qoo in the (Pc, Q) plan. The legend 
indicates the value of Qoo. The scalar held is no longer real 
for Q > l/(127r). All curves, except for Qoo = 0, reach a 
maximal value. 


We note that where the field becomes imaginary, its 
kinetic term has the wrong sign in the action. In GR 
this would inevitably lead to classical and quantum in¬ 
stabilities. However, in the present case, the situation 
is less clear due to the coupling of the kinetic term to 
the Einstein tensor. The stability of these configurations 
with imaginary field is an interesting and open question. 

Conclusion. — In this paper, we investigated the ef¬ 
fect of a scalar field with a non-minimal kinetic coupling 
to gravity. This model is a sector of the Horndeski’s the¬ 
ory of gravity and has a deep connection with Galileon 
gravity. We focused on the case with vanishing cosmo¬ 
logical constant and vanishing standard kinetic term cou¬ 
pling, since it admits exact vacuum Schwarzschild solu¬ 
tions. When the scalar field is linear in time we find that 
neutron stars exist for some range of the only free param- 


FIG. 4. Plot of F'^ for rj > 0 (left) and rj < 0 (right). In 
the case ri < 0, F' is imaginary close to the origin for small 
values of Q^, and becomes real for larger values of r, while 
it is always real for > 0. 


eter of the theory The equation of state used here 

captures very well the main effects of the derivative cou¬ 
pling. In particular, the reality conditions on the scalar 
field are independent of the equation of state so we do 
not expect large deviations if one chooses more realistic 
ones. 


We derived a constraint on r/Q^ from the existence of 
compact configurations, and argued that less dense stars, 
such as white dwarfs, are supported by this model. To 
date, this is the first constraint and astrophysical via¬ 
bility check of this model. The parameter Q is clearly 
non-perturbative since the “extremal” case Q = 0 still 
shows a deviation from GR. In general, when Q 0 the 
scalar field grows linearly in time but, since the dynamics 
depends only on its gradient, there are no diverging phys¬ 
ical observables. In fact, the scalar field plays the mere 
role of a “clock”, just like the scale factor in cosmology. 

A very attractive feature of our model is that the exter¬ 
nal metric is the same as in GR. From an observational 
point of view, this means that Horndeski’s stars pass all 
Solar System tests. In this respect, the model is similar 
to Palatini /(i?) or Eddington-inspired Born Infeld mod¬ 
els j23 (with similar effects on compact configurations 
[ 29 I l3Cll|). In fact, these theories are also indistinguish¬ 
able from GR in vacuum, but they are modified in the 
presence of matter fields 311. 

In summary, we show that ours is a sound astrophysi¬ 
cal model, which requires a minimal modification of GR 
and that does not violate any Solar System test. The un¬ 
derlying Horndeski’s gravity is currently under thorough 
scrutiny as a possible candidate to explain fundamental 
problems, such as inflation and dark energy. Therefore, 
we think that also the Horndeski’s star deserves full at¬ 
tention as it can shed some light on the general theory. 
We think that, to further asses the viability of our model, 
it would be very interesting to address questions such as 
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the the gravitational wave spectrum during a merger of 
Horndeski’s neutron stars or the observational signature 
on the pulsar emission in the presence of strong magnetic 
fields, and we hope to report soon on these issues. 
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